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ABSTRACT 

We provide X-ray constraints and perform the first X-ray spectral analyses for bright SCUBA sources 
(f&50pm > 5 mJy; S/N> 4) in an 8^4 x 8?4 area of the 2 Ms Chandra Deep Field-North survey containing the 
Hubble Deep Field-North. X-ray emission is detected from seven of the ten bright submm sources in this region 
down to 0.5-8.0 keV fluxes of « 1 x 10 1 6 erg cm 2 s 1 , corresponding to an X-ray detected submm source 
density of 36OIJ3Q deg 2 ; our analyses suggest that this equates to an X-ray detected fraction of the bright submm 
source population of £ 36%, although systematic effects may be present. Two of the X-ray detected sources 
have nearby (within 3") X-ray companions, suggesting merging/interacting sources or gravitational lensing ef¬ 
fects, and three of the X-ray detected sources lie within the approximate extent of the proto-cluster candidate 
CXOHDFN J123620.0+621554. Five of the X-ray detected sources have flat effective X-ray spectral slopes 
(r < 1.0), suggesting obscured AGN activity. X-ray spectral analyses suggest that one of these AGNs may be 
a Compton-thick source; of the other four AGNs, three appear to be Compton-thin sources and one has poor con¬ 
straints. The rest-frame unabsorbed X-ray luminosities of these AGNs are more consistent with those of Seyfert 
galaxies than QSOs (i.e.. Lx ~ 10 43 —10 44 erg s -1 ). Thus, the low X-ray detection rate of bright submm sources by 
moderately deep X-ray surveys appears to be due to the relatively low luminosities of the AGNs in these sources 
rather than Compton-thick absorption. A comparison of these sources to the well-studied heavily obscured AGN 
NGC 6240 shows that the average AGN contribution is negligible at submm wavelengths. The X-ray properties of 
the other two X-ray detected sources are consistent with those expected from luminous star formation; however, 
we cannot rule out the possibility that low-luminosity AGNs are present. The three X-ray undetected sources 
appear to lie at high redshift (z > 4) and could be either AGNs or starbust galaxies. 

Subject headings: cosmology: observations - galaxies: active - submillimeter - surveys - X-rays 


1. INTRODUCTION 

Due to the strong negative /("-correction for galaxies at sub¬ 
millimeter (submm; X =300-1000 fi m) wavelengths, submm 
surveys facilitate observations of the high-redshift universe 
(e.g., Blain & Longair 1993). For typical starburst galaxies and 
Active Galactic Nuclei (AGNs), the expected submm flux den¬ 
sity is approximately insensitive to redshift for z ~ 1-10. In¬ 
deed, the majority of the submm-detected sources in SCUBA 
(Holland etal. 1999) surveys appear to be extremely luminous 
infrared galaxies at z « 1-6 (e.g., Smail, Ivison, & Blain 1997; 
Hughes etal. 1998; Ivison etal. 1998; Smail etal.2000, 2002; 
Dunlop 2001; Blain etal. 2002). 

It is generally believed that these submm-detected sources 
are proto-galaxies undergoing intense dust-enshrouded star for¬ 
mation (e.g., Blain etal. 1999; Lilly etal. 1999). However, 
evidence for AGNs has also been found using a variety of 
techniques. For example, optical spectroscopic observations 
have revealed AGNs in a number of optically bright submm 
sources (e.g., Ivison etal. 1998; Barger etal. 1999a; Soucail 
et al. 1999), and SCUBA observations of known optically bright 
Quasi-Stellar Objects (QSOs; i.e., luminous type 1 AGN) have 
shown that luminous AGNs comprise a fraction of the submm 
source population (e.g., McMahon etal. 1999; Page etal. 2001; 
Isaak etal. 2002; Willott etal. 2002). 1 Although clearly some 
SCUBA sources host AGNs, the ubiquity of AGNs in the 
submm source population, and their contribution to the submm 


emission, are poorly constrained. 

Arguably the best discriminator of AGN activity is the de¬ 
tection of hard X-ray emission (i.e., > 2 keV). Hard X-ray 
observations are particularly useful in identifying AGNs in 
sources where the optical signatures are weak (e.g., optically 
faint AGNs or AGNs without strong emission lines; Vignati 
et al. 1999; Mushotzky etal. 2000; Alexander et al. 2001; Horn¬ 
schemeier etal.2001). Due to their high-energy X-ray cover¬ 
age, high X-ray sensitivity, and excellent positional accuracy, 
the Chandra (Weisskopf et al. 2000) and XMM-Newton (Jansen 
etal.2001) observatories offer the best opportunities for the 
X-ray investigation of SCUBA sources. Although sensitive 
enough to detect luminous AGNs out to high redshift, the cross¬ 
correlation of moderately deep X-ray surveys with SCUBA sur¬ 
veys has yielded little overlap between the X-ray and submm 
detected source populations (Fabian et al. 2000; Hornschemeier 
et al. 2000,2001; Severgnini et al. 2000; Barger et al. 2001a; Al- 
maini et al. 2002; however, see also Bautz et al. 2000 and Ivison 
et al. 2002). These studies have suggested that a large number 
of bolometrically important AGNs (i.e., AGNs that can dom¬ 
inate the submm emission) can only be present in the submm 
source population if they are obscured by Compton-thick (i.e., 
Ah > 1.5 x 10 24 cm 4 ) material. 

Targeting X-ray sources detected in the 1 Ms Chandra 
Deep Field-North (CDF-N; Brandt etal.2001) survey with 
SCUBA jiggle-map observations, Barger et al. (2001b) showed 
that submm counterparts can be found for faint X-ray sources. 


1 The exact fraction of optically bright QSOs in the submm source population is poorly known. Since only one of the 15 sources in the SCUBA-lens survey (Smail 
etal. 2002) is a QSO [SMM J02399-0136 is a broad absorption line QSO (BALQSO); see Vemet & Cimatti 2001], it is probably of the order of « 5-10%. 
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Fig. 1.— Adaptively smoothed full-band 2 Ms Chandra image. This image has been made using the standard ASCA grade set and has been adaptively smoothed 
at the 2.5o level using the code of Ebeling, White, & Rangarajan (2002). The HDF-N is shown as the polygon at the center of the image, and the large box indicates 
the reduced Hawaii flanking-field area used in this study. The circles indicate the source positions of the labelled submm sources (see §2.2 and Table 1). The dashed 
box indicates the position and approximate extent of CXOHDFN J 123620.0+621554 (e.g., Bauer et al. 2002a), a candidate high-redshift cluster or proto-cluster that 
has an overdensity of submm sources (see §2.4.1); since this source emits predominantly in the soft band it is not apparent in this image. Most of the faint extended 
emission is instrumental background or artifacts of the smoothing algorithm. The two white streaks towards the top of the image are due to gaps in the ACIS-I 
CCDs. The faintest X-ray sources are below the significance level of the smoothing and thus are not visible in this figure. 


Five of the six X-ray detected submm sources with good X-ray on the origin of the submm emission from AGN-classified 

constraints had flat (r < 1.0) effective X-ray spectral slopes, submm sources by determining if their AGNs are obscured 

an indicator of obscured AGN activity. However, while clearly by Compton-thin or Compton-thick material. In this study 

showing that many of the AGNs found in submm sources are we focus on an 84 x 84 region within the 2 Ms CDF-N sur- 

obscured, Barger et al. (2001b) did not investigate the possibil- vey with excellent multi-wavelength data, including deep radio 

ity that the absorption toward the AGNs is Compton thick. Thus (Richards et al. 1998; Richards 2000) and optical/near-infrared 

strong constraints on the origin of the submm emission were not (near-lR) observations (Barger etal. 1999b). Due to the ex- 

placed. tremely high sensitivity of the 2 Ms CDF-N survey, we have 

The primary aim of this paper is to place strong constraints the potential to detect X-ray emission from high-luminosity 
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starburst galaxies (i.e.. Lx ~ 10 43 erg s" 1 ) to z ~ 2.5. In this 
study we place the tightest constraints to date on the X-ray 
emission from star formation in submm-detected galaxies. The 
Galactic column density along this line of sight is (1 .6 ± 0.4) x 
10 20 cm -2 (Stark etal. 1992), and Hq = 65 km s -1 Mpc~\ 
Bm = 5 , and O a = | are adopted throughout. 

2. OBSERVATIONS AND SOURCE SAMPLES 

2.1. X-ray Observations 

The X-ray results reported in this paper were obtained 
with ACIS-I (the imaging array of the Advanced CCD Imag¬ 
ing Spectrometer; Garmire et al. 2002) onboard Chandra. 
The 2 Ms CDF-N observations were centered on the Hub¬ 
ble Deep Field-North (HDF-N; Williams etal. 1996) and cover 
« 460 arcmin 2 . The region investigated in this study extends 
over an 8?4 x 8 '.4 area that has its center near to the aim-point 
of the CDF-N observations. This 70.3 arcmin 2 * region is defined 
by the optical and near-IR observations of Barger et al. (1999b) 
and is referred to here as the reduced Hawaii flanking-field area 
(see Alexander et al. 2001). The X-ray data processing of these 
observations was similar to that described in Brandt et al. (2001) 
for the 1 Ms Chandra exposure. The full 2 Ms source catalog 
is presented in Alexander et al. (2003). 

One hundred and ninety-three (193) X-ray sources are de¬ 
tected in the reduced Hawaii flanking-field area with a WAVDE- 
TECT (Freeman et al. 2002) false-positive probability threshold 
of 10 - 7 down to on-axis 0.5-2.0 keV (soft-band) and 2-8 keV 
(hard-band) flux limits of « 1.5 x 10 ~ 17 erg cm ~ 2 s -1 and ss 
1.0 x 10 16 erg cm ~ 2 s ' 1 , respectively. The X-ray coverage in 
this region is fairly uniform, with 0.5-8.0 keV (full-band) effec¬ 
tive exposure times of the 193 X-ray detected sources covering 
1.28-1.94 Ms (with a median time of 1.82 Ms). The adaptively 
smoothed full-band 2 Ms Chandra image is shown in Figure 1. 
The positional uncertainties of the X-ray sources are 4S 1 "0, 
depending on off-axis angle (see Alexander et al. 2003), and the 
median positional uncertainty is 0."4. These small positional 
uncertainties allow for accurate cross-correlation with multi¬ 
wavelength counterparts. 

2.2. Bright Submm Sources 

In this study we focus on bright submm (fysopm > 5 mJy) 
sources detected with a signal-to-noise ratio (S/N) > 4. 
This threshold was chosen to guard against possible spurious 
SCUBA sources (e.g., see the estimated number of spurious 
sources at different significance thresholds for the SCUBA ob¬ 
servations of the ELAIS-N2 and Lockman Hole East regions; 
Scott et al. 2002). 

The submm data used here was taken from a number of dif¬ 
ferent studies (Hughes et al. 1998; Barger et al. 2000, 2002; Bo- 
rys etal.2002). Ten bright submm sources are detected with 
S/N> 4; see Table 1 and Figure 1 ? If we relaxed the detection 
threshold to S/N> 3.5 another five sources would be added to 
our sample; however, « 20 % of the sources in the total sam¬ 
ple could be spurious (e.g., Scott etal.2002). Although we 
have utilized all of the published SCUBA data currently avail¬ 
able, this region only has complete SCUBA coverage down 
to fssOfjm = 12 mJy (« 75% of the region has SCUBA cov¬ 
erage down to fssOfjm = 5 mJy; e.g., see Figure 1 of Barger 


et al. 2002). Therefore, there are likely to be further bright 
submm sources in this region that have not yet been detected 
by SCUBA. 

Two main observing techniques were employed in these 
studies: Hughes etal. (1998) and Borys etal. (2002) performed 
SCUBA scan-map observations while Barger etal. (2000, 
2002) performed SCUBA jiggle-map observations of radio and 
X-ray sources. The positional uncertainty of a SCUBA source 
is dependent on the signal-to-noise ratio and is typically « 2-A" 
(e.g., Serjeant et al. 2002; Smail et al. 2002; although see Hogg 
2001). However, in the cases where radio and X-ray sources 
have been targeted with SCUBA jiggle-map observations, <\" 
positions can be inferred from the radio and X-ray data. 

Seven of the ten submm sources were detected in SCUBA 
jiggle-map observations of radio and X-ray sources, and have 
positional uncertainties of <\". For these sources we identified 
optical and near-IR counterparts in the photometric source cat¬ 
alogs produced by Alexander et al. (2001) using a matching 
radius (see Table 1); these source catalogs were produced from 
the Barger etal. (1999b) /-band and /L/C-band images which 
reach « 2o magnitude limits of I = 25.3 and IIK' = 21.4, re¬ 
spectively . 3 

Three of the ten submm sources were detected in SCUBA 
scan maps and have positional uncertainties of « 2^1". How¬ 
ever, in the case of one source (123652.0+621225, which is of¬ 
ten referred to as HDF850.1, the brightest submm source in the 
HDF-N; Hughes etal. 1998) we can adopt the sub-arcsecond 
1.3 mm source position obtained with the IRAM Interferome¬ 
ter on the Plateau de Bure (IRAM PdBI) telescope (i.e., Downes 
etal. 1999). Dunlop etal. (2002) used these data in combi¬ 
nation with extremely sensitive MERLIN+VLA 1.4 GHz ra¬ 
dio observations to identify accurately the optical and near- 
IR counterpart for 123652.0+621225. Since this source is lo¬ 
cated in the HDF-N, the optical and near-IR constraints are 
more sensitive than for the other submm sources investigated 
here (123652.0+621225 has I > 28.6 and K = 23.4; Dunlop 
et al. 2002). We must rely on the SCUBA positions for the other 
two sources since neither source has a radio or X-ray counter¬ 
part within 4". The « 2—4" positional uncertainties for these 
sources are too large to be able to identify optical and near-IR 
counterparts relibably. 

2.3. Source Redshifts 

Only one of the submm sources has an optical spectroscopic 
redshift (123629.1+621045 has z = 1.013; Cohen etal.2000; 
Hornschemeier etal.2001) and another source has a redshift 
estimate based on its multi-wavelength spectral energy dis¬ 
tribution (SED; 123652.0+621225 has z = 4.1+^; Dunlop 
et al. 2002). However, on the assumption that the other sources 
obey the radio-to-far-infrared (far-IR; A. =40-120 qtu) correla¬ 
tion found for local star-forming galaxies (e.g., Helou, Soifer, & 
Rowan-Robinson 1985), we can estimate redshifts for all of the 
radio-detected sources, and place redshift lower limits on the 
radio-undetected sources, using the millimetric redshift tech¬ 
nique (e.g., Carilli & Yun 1999). Since the submm emission 
corresponds to the longer wavelength Rayleigh-Jeans tail of the 
thermal dust emission detected at far-IR wavelengths, redshifts 
are directly determined (albeit with considerable uncertainty) 
on the basis of the radio-to-submm spectral slope. As noted by 


2 We also included the studies of Chapman et al. (2001) and Serjeant et al. (2002) but no further SCUBA sources met our criteria. 

3 These images are publicly available at http://www.ifa.hawaii.edu/~cowie/hdflank/hdflank.html. The relationship between the K band and HK' band is HK' — K = 

0.13 + 0.05(7 —AT) (Barger etal. 1999b). 



4 


ALEXANDER ET AL. 


Smail et al. (2000), this technique should give a lower limit on 
the source redshift when a component of the radio emission is 
produced by an AGN. 

We determined millimetric redshifts using the redshift es¬ 
timators of Carilli & Yun (2000), Eales et al. (2000), and 
Rengarajan & Takeuchi (2001). We found that the Eales 
et al. ( 2000 ) estimator gave the lowest redshift determinations 
(a minimum redshift of z = 1.3 and a maximum redshift of 
z > 3.2), while the Rengarajan & Takeuchi (2001) estimator 
gave the highest redshift determinations (a minimum redshift 
of z = 2.0 and a maximum redshift of z > 4.8). The differences 
in the redshift determinations are due to the assumed SED and 
the temperature of the dust that produces the far-IR-to-submm 
emission. The millimetric redshifts of the source with an opti¬ 
cal spectroscopic redshift (123629.l+621045;z = 1.013) cover 
Z = 1.4-2.3 (z = 1.0-3.0 when the la millimetric redshift un¬ 
certainties are taken into account). Clearly the optical spec¬ 
troscopic redshift of 123629.1+621045 is marginally in agree¬ 
ment with the millimetric redshift estimates (we note that larger 
studies find a ^ 50% agreement between the optical spectro¬ 
scopic redshifts and millimetric redshifts of submm sources; 
e.g., Smail etal. 2000, 2002). However, in the absence of better 
redshift estimations we must determine the redshifts for the ma¬ 
jority of these sources using the millimetric redshift technique. 

The adopted millimetric redshifts for our sources are deter¬ 
mined using the Carilli & Yun (2000) redshift estimator, which 
provides a good compromise between the Eales et al. (2000) 
and Rengarajan & Takeuchi (2001) millimetric redshift esti¬ 
mates (see Table 1). We note that since the typical la millimet¬ 
ric redshift uncertainty for a source is A z = 0.8 (the range of 
redshift uncertainties is A z = 0 .3-1.1), the differences in these 
redshift determinations are not critical. However, even given 
these uncertainties, all sources have z > 1 with all millimetric 
redshift estimators. 

2.4. X-ray Detected Submm Sources 

Five of the ten submm sources were detected by Barger 
etal. (2002) in SCUBA jiggle-map observations of X-ray 
sources and are thus X-ray detected submm sources (see Ta¬ 
ble 2). We searched for X-ray counterparts for the other five 
submm sources using a matching radius for the three sources 
with a radio counterpart and a A" matching radius for the two 
sources without a radio counterpart . 4 One further X-ray coun¬ 
terpart was found for one of the radio-detected submm sources 
(123622.6+621629; see Table 2). 

Given the low surface density of the submm sources, we 
can also search for lower significance X-ray sources associated 
with submm sources without introducing a significant number 
of spurious X-ray sources. We ran WAVDETECT with a false¬ 
positive probability threshold of 10 5 and found one further 
match to a radio-detected submm source (123618.3+621551; 
see Table 2); the probability of this X-ray-submm match being 
spurious is < 1%. This source was also detected in the 1 Ms 
Chandra exposure at a 10 5 false-positive probability thresh¬ 
old (Alexander etal.2001) and is probably not detected here 
at a higher probability threshold since it lies within the bright 
extended X-ray source CXOHDFN J 123620.0+621554 (e.g., 
Bauer etal. 2002a). Thus, seven (70^%) of the ten submm 
sources have X-ray counterparts . 5 

Three submm sources are X-ray undetected. There is no evi¬ 


dence for excess X-ray emission in the Chandra images for two 
of these sources (123620.3+621701 and 123621.1+621250); 
however, there is a suggestion of soft-band emission at the 
location of 123652.0+621225. Indeed, 123652.0+621225 lies 
within 0"2 of the position of an X-ray source detected by 
WAVDETECT at a 10 4 false-positive probability threshold in 
the 1 Ms CDF-N study of Very Red Objects (VROs, I — K> 4; 
Alexander etal.2002b). Since 123652.0+621225 lies closer 
to the X-ray source than the VRO reported in Alexander 
et al. (2002b), it may actually be the correct counterpart to this 
extremely faint low-significance X-ray source. However, since 
the X-ray source is not detected at a higher significance level 
in the 2 Ms Chandra exposure, we consider 123652.0+621225 
undetected here. 

Examining the association between Chandra and SCUBA 
sources, Almaini et al. (2002) found little overlap between the 
source populations but detected a strong clustering signal. They 
suggested that Chandra and SCUBA sources are generally un¬ 
related (at least down to the depths of their 75 ks Chandra sur¬ 
vey) but trace the same large-scale structure. Although there 
may be some large-scale structure effects present in our sam¬ 
ple, the cross-correlation distances between our Chandra and 
SCUBA sources are much smaller than the up to 100 " cluster¬ 
ing distances found by Almaini etal. (2002). Thus, we expect 
all of our X-ray-submm matches to be legitimate. 

2.4.1. X-ray Detected Submm Source Density 

The derived density of X-ray detected submm sources in this 
region with /sso^m > 5 mJy is 360^3^ deg -2 . This source den¬ 
sity should be considered a lower limit since there may be fur¬ 
ther bright submm sources with X-ray counterparts in this re¬ 
gion that have not yet been detected by SCUBA (see §2.2). 

We can compare our source density to that found in other 
studies. The SCUBA survey of the ELAIS-N2 and Lock- 
man Hole East regions (i.e., Scott etal.2002) contains sensi¬ 
tive Chandra and XMM-Newton observations (down to 0.5- 
8.0 keV flux limits of « 10 -15 erg cm -2 s -1 ; i.e.. Manners 
et al. 2002; Hasinger et al. 2001). One Chandra counterpart was 
matched to a submm source in the ELAIS-N2 region (Almaini 
et al. 2002), corresponding to an X-ray detected submm source 
density of 35^ deg 2 . Four XMM-Newton counterparts were 
found for submm sources in the Lockman Hole East region 
(Ivison et al. 2002), corresponding to a higher X-ray detected 
submm source density of 120+^ deg -2 . The X-ray detected 
submm source densities found in our study are « 3-10 times 
higher than those found by Almaini etal. (2002) and Ivison 
et al. (2002) due to the greater sensitivity of the 2 Ms Chandra 
observations. 

Three X-ray detected submm sources lie within 1' of the X- 
ray source CXOHDFN J123620.0+621554 (see Figure 1). This 
extended X-ray source could be a high-redshift (z ~ 1-2) clus¬ 
ter or proto-cluster (Bauer et al. 2002a; Chapman et al. 2002b), 
and the similar millimetric redshifts of the three nearby X-ray 
detected submm sources suggest that they may be directly asso¬ 
ciated with CXOHDFN J 123620.0+621554. The statistics are 
extremely poor on how common such sources are. However, if 
they are rare, then our X-ray detected submm source density is 
likely to be higher than average at this X-ray depth; thus our 
sample may be affected by “cosmic variance”. 


X-ray sources detected in the soft band (0.5-2.0 keV), hard band (2-8 keV) or full band (0.5-8.0 keV) were matched to submm counterparts. 
All errors are taken from Tables 1 and 2 of Gehrels (1986) and correspond to the lo level; these were calculated assuming Poisson statistics. 
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Fig. 2.— /-band images of the X-ray detected submm sources with nearby X-ray companions. The contours indicate the extent and intensity of the X-ray 
emission; the lowest contour corresponds to 2 counts pixel -1 and each successive contour represents a factor of two increase in the number of counts pixel -1 . The 
large crosses mark the locations of the 1.4 GHz radio counterparts and the assumed submm counterparts. Each image is 9”43 on a side. 


2.4.2. Fraction of Submm Sources with X-ray Counterparts 

The reduced Hawaii flanking-field area does not have com¬ 
plete SCUBA coverage down to /g 50 /jm = 5 mJy (see §2.2), and 
hence we cannot directly place constraints on the fraction of 
the bright submm source population with X-ray counterparts. 
However, based on the submm source density results of Eales 
et al. (2000), and the compilations of Blain et al. (2002), Cowie, 
Barger, & Kneib (2002), and Smail etal. (2002), we would ex¬ 
pect 300-1000 submm sources per deg 2 with fssQ^m > 5 mJy. 
This suggests that the fraction of bright submm sources in this 
region with X-ray counterparts is likely to be x) 36%. Note that 
this source fraction was not directly measured in our field and 
systematic effects could be present. 

2.4.3. Nearby X-ray Companions 

Two of the X-ray detected submm sources 
(123616.1+621513 and 123646.0+621448) have X-ray coun¬ 
terparts within 3" (see Figure 2; Alexander et al. 2003). In both 
cases one X-ray source has a radio counterpart and one does 
not. In addition, the radio-detected sources have I « 24-25 
while the radio-undetected sources are optically blank. Given 
the density of X-ray sources in this region, the probability of 
two X-ray sources lying within 3" of each other is f 1 % and 
therefore these associations appear to be significant. 

There are two obvious scenarios that can explain the close 
proximity of the sources in these two pairs: (1) these sources are 
merging or interacting, or (2) these sources are gravitationally 
lensed. Unfortunately, since these sources are optically faint, 
we cannot search for the features of merging, interacting, or 
gravitational lensing in the current optical images. However, 


in the gravitational-lensing scenario, each member of each pair 
would generally have similar fluxes at all wavelengths (e.g., see 
review by Blandford & Narayan 1992). Therefore, the fact that 
one member of each pair is detected at both optical and radio 
wavelengths, while the other is not, provides some evidence 
against gravitational lensing. 6 

In both cases we have identified the radio-detected X-ray 
source as the counterpart to the submm source since the radio 
and submm emission are more closely matched in frequency 
space than the X-ray and submm emission. However, it is pos¬ 
sible that the radio-undetected X-ray sources could contribute 
to a fraction of the submm emission. 


3. X-RAY PROPERTIES OF SUBMM-DETECTED SOURCES 
3.1. Basic Source Properties 

We show the submm flux density versus full-band X-ray flux 
of all the sources in Figure 3. As a comparison we also show the 
submm-detected ROSAT-selected optically bright QSOs (i.e., 
luminous type 1 AGNs) of Page etal. (2001). There is a clear 
distinction in X-ray flux between the QSOs and our X-ray de¬ 
tected submm sources. The former would be detectable in 
;$ 10 ks Chandra exposures, while the majority of the latter are 
only detectable in f 100 ks Chandra exposures. Since optical 
spectroscopic observations of submm sources have shown that 
optically bright QSOs probably comprise only a small fraction 
of the submm source population (i.e., « 5-10%, see footnote 1; 
e.g., Barger etal. 1999a; Ivison etal.2000; Smail etal.2002), 
Ai 100 ks Chandra exposures are clearly required to detect X- 
ray counterparts for typical bright submm sources. 


6 The deep and high-spatial resolution Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations to be taken in this field as part of the Great 
Observatories Origins Deep Survey (GOODS) project should provide a direct test of these two scenarios. See http://www.stsci.edu/science/goods/ for details of the 
GOODS project. 
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Fig. 3.— Submm flux density versus full-band X-ray flux. The filled cir¬ 
cles are the X-ray detected submm sources, and the open circles are the X-ray 
undetected submm sources (see §2.4); crosses indicate sources classified as 
AGNs, and “U” indicates sources that do not have good X-ray spectral con¬ 
straints (see §3.2 and Table 2). The filled triangles are the submm-detected 
QSOs of Page etal. (2001); see §3.1. The upper limit arrows in the y-axis 
direction show the approximate full-band flux limits for Chandra ACIS-I sur¬ 
veys of different exposure times. While submm-detected QSOs are detectable 
in £ 10 ks Chandra exposures, ^ 100 ks Chandra exposures are required to 
detect typical submm sources. 

Figure 3 shows evidence for an anti-correlation between the 
full-band flux and the submm flux density for our sources. 
While this anti-correlation could be real, we caution that it may 
be an artifact caused by small number statistics and sample in¬ 
completeness. 

3.2. X-ray Source Classification 
3.2.1. X-ray Spectral Slopes and Luminosities 

The effective X-ray photon index (T) of a source can provide 
a simple constraint on its nature . 7 Within the context of AGN- 
classified sources, obscured AGNs have considerably flatter ef¬ 
fective X-ray spectral slopes than the canonical T ft 2.0 photon 
index of unobscured AGNs (e.g., George et al. 2000), due to the 
energy-dependent photo-electric absorption of the X-ray emis¬ 
sion (e.g., Risaliti, Maiolino, & Salvati 1999). Star-forming 
galaxies are also distinguishable from obscured AGNs since 
their rest frame ft 2-8 keV emission is often consistent with 
that of a r « 2.0 power law (e.g., Kim, Fabbiano, & Trinchieri 
1992a; Ptak et al. 1999). 8 

The X-ray luminosity of a source can provide a fur¬ 
ther check on the nature of the X-ray emission since very 
few starburst galaxies have full-band luminosities in ex¬ 
cess of Lx ~ 10 42 erg s -1 , even when including luminous 
sources at moderate redshifts (e.g., Moran, Lehnert, & Helfand 


Fig. 4.—X-ray band ratio versus redshift. The symbols have the same mean¬ 
ing as in Figure 3; the X-ray detected submm source 123618.3+621551 cannot 
be plotted here since it is only detected in the full band. The error bars in the 
jt-axis direction correspond to the uncertainties in the millimetric redshifts, and 
the error bars in the y-axis direction correspond to the uncertainties in the band 
ratio. The dashed lines show the expected band ratios for different absorb¬ 
ing column densities with an underlying Y = 2.0 power-law continuum; these 
have been calculated using PIMMS Version 3.2d. The AGN-classified sources 
are obscured, with implied intrinsic column densities of /Vj[ > 10 23 cm -2 . 

1999; Zezas, Alonso-Herrero, & Ward 2001; Alexander 
etal.2002a). While this would usually be a reasonable up¬ 
per limit to the expected X-ray emission from star-formation 
activity, given the extreme infrared luminosities of submm 
sources (Lir = 10 12 -10 13 L 0 ; e.g., Ivison etal. 1998; Smail 
et al. 2002 ), we might expect full-band luminosities up to an or¬ 
der of magnitude larger (i.e., up to Lx ft 10 43 erg s' 1 ) for this 
source class . 9 

The X-ray band ratio versus redshift for the six sources with 
band ratio constraints is shown in Figure 4. Five of the X-ray 
detected submm sources have effective photon indices of T < 
1.0. Such flat X-ray spectral slopes are only likely to be pro¬ 
duced by obscured AGN activity and, under the assumption that 
the underlying continuum follows a T = 2.0 power law, these 
sources have intrinsic column densities of Ah ^ 10 23 cm“ 2 . 

The rest-frame X-ray luminosities of these sources are calcu¬ 
lated as 

L x = 4ndlf x (1 +z) r “ 2 erg s _1 , (1) 

where d l is the luminosity distance (cm), fx is the X-ray flux 
(erg cm -2 s -1 ), and T = 2.0 is the assumed photon index. 

Assuming the underlying continuum follows a T = 2.0 
power law, the unabsorbed rest-frame full-band luminosities 
(Lx ~ 10 43 —10 44 erg s -1 ) of the five sources with flat X-ray 


7 The photon index is related to the energy index by T = a+1 where F v v a . 

8 Detailed X-ray studies of local star-forming galaxies have shown that they have more complex X-ray spectra than simple power-law emission (e.g., power-law 
emission and a very soft ~ 0.7 keV thermal component; Ptak etal. 1999); however, at the probable redshifts of our sources the very soft thermal component will be 
redshifted out of the Chandra energy band. 

9 While Lir = 10 12 -10 13 L© galaxies are probably the best candidates for detecting Lx ~ 10 43 erg s -1 starburst galaxies, none have been found to date. 
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spectral slopes are also consistent with obscured AGN activity 
(the intrinsic absorption is calculated from the band ratios; i.e., 
see Figure 4). These sources are classified hereafter as AGNs 
(see Table 2). 

The X-ray spectral constraints are weak for the other five 
submm sources. The rest-frame full-band luminosities or up¬ 
per limits (Lx ~ 10 42 -10 43 erg s _1 ) for the three sources with 
luminosity constraints (see Table 2) suggest that their X-ray 
emission could be produced by luminous star-formation activ¬ 
ity. However, since we cannot rule out AGNs in any of these 
five sources (e.g., they could be extremely highly obscured or 
of low luminosity), they are classified here as unknown (see 
Table 2). 

3.2.2. X-ray Emission from Star Formation 

From a comparison of the X-ray and radio luminosities of 
starburst galaxies detected in the 1 Ms Chandra exposure of 
the CDF-N, Bauer et al. (2002b) showed that the radio emission 
can be used to predict the X-ray emission from star formation 
(see also Ranalli, Comastri, & Setti 2002). The main assump¬ 
tion in this prediction is that the radio emission is dominated by 
star formation; the presence of an AGN component to the radio 
emission will lead to an overprediction of the contribution from 
star formation at X-ray energies. 

None of the radio-detected submm sources is classified by 
Richards (1999) as an AGN on the basis of its radio emission 
although four sources are classified as unknown (see Table 1). 
Since only two of these four sources are classified as AGNs at 
X-ray energies, and all but one have steep radio spectral slopes, 
the radio emission is likely to be dominated by star-formation 
activity in most cases (e.g., Condon 1992; Richards 2000). The 
rest-frame radio luminosity densities of these sources are cal¬ 
culated as 

Lx A GHZ = 4 7t d\ f\ 4 GHz 10^ (1+z) a “' W Hz^, (2) 

where d\ is the luminosity distance (cm), / 1.4 GHz is the 
1.4 GHz flux density (/By), and a corresponds to the radio spec¬ 
tral index. In determining the radio luminosity density we have 
assumed a = 0.8, the average spectral index for star-forming 
galaxies (e.g., Yun, Reddy, & Condon 2001). 

A comparison of the rest-frame full-band X-ray and radio 
luminosities of the radio-detected submm sources is shown in 
Figure 5. The predicted X-ray luminosity from star-formation 
activity for all of the radio-detected submm sources is Lx ~ 
10 42 -10 43 erg s" 1 . The AGN-classified sources are on average 
« 10 times more luminous at X-ray energies than that predicted 
by star formation. However, the X-ray emission from the other 
three radio-detected submm sources could be completely dom¬ 
inated by star-formation activity. 

3.3. X-ray Spectral Analysis 

Our interpretation of the X-ray spectral slopes and luminosi¬ 
ties for the five AGN-classified sources (see §3.2.1) was based 
on the assumption that the obscuration to the X-ray emitting 
source is not optically thick to Compton scattering. How¬ 
ever, in some cases the obscuration may be Compton thick 
(i.e., Ah > 1.5 x 10 24 cm -2 ), and the observed X-ray emis¬ 
sion will be predominantly produced by reflection and scatter¬ 
ing processes (e.g.. Matt etal. 1997; Turner etal. 1997a; Matt 



Fig. 5.— Rest-frame full-band X-ray luminosity versus rest-frame 1.4 GHz 
radio luminosity density for the radio-detected submm sources. The radio- 
undetected sources cannot be plotted as they all have redshift lower limits. 
The symbols have the same meaning as in Figure 3. The error bars correspond 
to the uncertainties in the luminosities. The X-ray luminosities of the AGN- 
classified sources correspond to the unabsorbed luminosity. The shaded region 
denotes the 1 o dispersion in the locally-determined X-ray/radio correlation 
(see Figure 6 of Shapley, Fabbiano, & Eskridge 2001). The dashed line is the 
regression fit to emission-line galaxies investigated in Bauer et al. (2002b) with 
the 1 Ms Chandra exposure of the CDF-N. The AGN-classified sources have 
X-ray luminosities on average « 10 times higher than those expected from star 
formation. The X-ray luminosities of the sources without good X-ray spectral 
constraints are consistent with those expected from star formation. 

et al. 2000). The most direct discrimination between Compton- 
thin and Compton-thick absorption is made with X-ray spectral 
analysis. The X-ray spectrum of a Compton-thick AGN is gen¬ 
erally characterized by a large equivalent width iron Ka emis¬ 
sion line (EW 0.5 keV; e.g., Matt, Brandt, & Fabian 1996; 
Maiolino etal. 1998; Bassani etal. 1999; Matt etal. 2000) and 
a flat or inverted (r < 0.0) X-ray spectral slope, due to pure 
reflection. By contrast, the X-ray spectrum of a Compton-thin 
AGN is usually well fitted by an absorbed power-law model 
and a smaller equivalent width iron Ka emission line (gener¬ 
ally EW ~ 0.1-0.5 keV; Nandraetal. 1997; Turner etal. 1997b; 
Bassani etal. 1999). 

To account for the range of roll angles and aim points in the 
20 separate observations that comprise the 2 Ms CDF-N, we 
used the ACIS source extraction code (ACIS Extract) de¬ 
scribed in Broos etal. (2002). 10 Briefly, for each source this 
code extracts the counts from each of the observations taking 
into account the changing shape and size of the PSF with off- 
axis angle as given in the Chandra X-ray Center (CXC) PSF 
library. 11 A local background is extracted after all sources are 
excluded from the X-ray event file, and the spectra and re¬ 
sponse matrices are summed using standard FTOOLS routines 
(Blackburn 1995). The latest version of XSPEC (vll.2.0i; Ar- 
naud 1996, 2002) was used for all of the model fitting, and 


10 ACIS EXTRACT is a part of the TARA software package and can be accessed from http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html. 

11 See http://asc.harvard.edU/ciao2.2/documents_dictionary.html#psf. 
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the fit parameter uncertainties are quoted at the 90% con¬ 
fidence level (i.e., Avni 1976). Given the limited counting 
statistics for the AGN-classified sources (« 30-200 full-band 
counts), we performed all the X-ray spectral analysis using 
the C-statistic (Cash 1979). The advantage of using the C- 
statistic is that the data can be fitted unbinned, making it well 
suited to low-count sources (e.g., Nousek & Shue 1989). In 
the latest version of XSPEC, the C-statistic can be used on 
background-subtracted data. We carried out several checks 
of the background-subtraction method to verify that no spuri¬ 
ous residual features were present in the background-subtracted 
data. 

We only considered two simple models for the X-ray spec¬ 
tral fitting: (1) a power-law model with Galactic absorption, and 
(2) a power-law with both Galactic and additional absorption. 
The results are given in Table 3. As expected, the photon in¬ 
dices derived using the Galactic absorbed power-law emission 
model are in good agreement (within ±0.2) with those calcu¬ 
lated from the band ratio. When fit with an absorbed power-law 
model, the best-fitting photon indices and observed absorption 
for three of the sources are in broad agreement with those ex¬ 
pected for Compton-thin AGNs (e.g.. Turner etal. 1997b; Bas- 
sani etal. 1999), while two sources have particularly flat pho¬ 
ton indices (123622.6+621629 and 123712.0+621211), pos¬ 
sibly indicative of reflection and Compton-thick absorption. 
However, with only « 30 full-band counts, 123712.0+621211 
is probably too faint to provide reliable constraints. 

To investigate whether any of these sources could be 
Compton-thick AGNs we searched for emission-line features in 
the 2 o binned X-ray spectra. The only source that showed ev¬ 
idence for an emission-line feature was 123622.6+621629 (see 
Figure 6); this was also one of the sources with a flat photon 
index. To assess whether this « 4.4 keV emission-line fea¬ 
ture was real, we simulated 10,000 low-count (« 90 full-band 
counts) X-ray spectra using the fakeit routine within XSPEC; 
the input model was a Galactic absorbed T = -0.7 power law. 
Thirty-four (0.34%) of the simulations produced a feature with 
an equivalent width as large as or larger than that measured in 
the data, corresponding to a significance level of 2.9 o. While 
this provides an estimation of the significance of this feature if 
we were expecting it to appear at« 4.4 ke V, it does not take into 
account of the fact that we simply identified the most significant 
feature over the observed energy range of « 1.0-6.4 keV. To 
take this into account we must determine the number of possible 
distinguishable features over this energy range. Given the aver¬ 
age spectral resolution of ACIS-I (« 150 eV) we would expect 
« 30 distinguishable features. However, by adding unresolved 
gaussian lines we estimate that we would only be able to distin¬ 
guish six features in the X-ray spectrum of 123622.6+621629; 
this demonstrates that photon statistics are more important than 
spectral resolution in low-count sources. Thus, the adjusted sig¬ 
nificance level for this feature is ss 2.3 o (i.e., 6 x 0.34%). 

3.3.1. 123622.6+621629: a Compton-thick AGN? 

Assuming that the « 2.3 o emission feature is a redshifted, 
unresolved 6.4 keV iron Ka emission line, the observed en¬ 
ergy (4.4±0.1 keV) implies z = 0.46 an d a rest-frame 

equivalent width of 0.71 qj keV. The moderately large equiv¬ 
alent width of this possible 6.4 keV iron Ka emission line, and 
the inverted X-ray spectral slope (r = -0.7; see Table 3), pro¬ 
vide evidence that 123622.6+621629 may be a Compton-thick 
AGN. Under the assumption that the underlying AGN model is 
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FIG. 6.— The X-ray spectrum of 123622.6+621629. The data are shown 
here in 2 o bins for presentation purposes; the error bars are Gaussian and un¬ 
derestimate the true errors by ~ 30%. The top panel shows the X-ray spectrum 
with a fitted model of Galactic absorbed power-law emission (T =-0.7; see Ta¬ 
ble 3). The bottom panel shows the ratio between the X-ray spectrum and the 
fitted model. Under the assumption that the « 2.3 o feature at 4.4 ±0.1 keV is 
the 6.4 keV iron Ka emission line, the inferred source redshift is z — 0.461q'q 2 
and the large rest-frame equivalent width (0.7 ^q ^ keV) provides evidence that 
123622.6+621629 may be a Compton-thick AGN (see §3.3.1). 

the same as that fitted to the Compton-thick AGN NGC 6240 
(Vignati etal. 1999), the unabsorbed rest-frame full-band lumi¬ 
nosity at z ~ 0.46 is Lx ~ 10 43 erg s -1 . 

The inferred redshift (z = 0 . 46 ^qq 2 ) is in considerable dis¬ 
agreement with the millimetric redshift (z = 2.41q g). However, 
we note that since 123622.6+621629 is detected in both the U 
band (Barger etal. 2000) and the U n band (Hogg etal.2000), 
it is actually likely to lie at z < 2 (to avoid the Lyman break). 
With a faint optical counterpart (I = 23.4), the absolute optical 
magnitude of 123622.6+621629 at z ~ 0.46 would correspond 
to that expected from a dwarf galaxy (i.e., w 3^1 mags below 
L , ). There are no known examples of luminous AGNs in dwarf 
galaxies in the local universe (e.g., McLeod & Rieke 1995), al¬ 
though it is possible that 123622.6+621629 is anL* galaxy with 
« 3-4 mags of dust extinction. 

If 123622.6+621629 lies at redshift of z « 0.46 then its prop¬ 
erties are atypical for a submm source (e.g., Ivison etal.2000; 
Smail etal.2002). However, we note that 123622.6+621629 
shares many similar characteristics with FN1-40, a z = 0.45 
submm source with an / = 24 optical counterpart (Chapman 
et al. 2002a). Although good X-ray constraints do not exist for 
FN1-40, the lack of broad optical emission lines led Chapman 
et al. (2002a) to claim no evidence for AGN activity. While 
FN1-40 is clearly not a broad-line AGN, we note that with an 
[OIII] A.5007/H(3 emission line ratio of 3 (see Figure 2 of 
Chapman et al. 2002a), it could be a narrow-line AGN (e.g., 
Baldwin, Phillips, & Terlevich 1981; Veilleux & Osterbrock 
1987). Hence, FN1-40 and 123622.6+621629 may be similar 
objects. Optical spectroscopy is clearly required to determine 
the redshift of 123622.6+621629 and help shed light on the ori¬ 
gin of its « 2.3 a emission line feature. 

4. THE ORIGIN OF THE SUBMM EMISSION 

It is often assumed that the submm emission from sources 
detected in SCUBA surveys can provide an unobscured indica- 
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Fig. 7.— Submm-to-X-ray spectral index (asx) versus redshift. The symbols have the same meaning as in Figure 3; the long-dashed line connecting two filled 
circles indicates the two possible redshifts for 123622.6+621629 (i.e., z = 0.46 1 ]]]]] and z = 2.4^,'g). The error bars in the x-axis direction correspond to the 
uncertainties in the millimetric redshifts, and the error bars in the y-axis direction correspond to the uncertainties in asx- The X-ray flux density corresponds to 
2 keV and is calculated from the full-band flux using the measured photon index of each source. Overplotted are the expected asx values for 3C273, NGC 6240 
(f sc = 0.05), and Arp 220 (dashed lines); alternative curves for NGC 6240 with less internal absorption ( Nh = 5 x I (A cm 2 ) and a smaller scattered flux fraction 
(fsc = 0.01) are also shown (dotted lines). This figure has been adapted from Fabian etal. (2000); the references for the 3C273, NGC 6240 and Arp 220 data are 
Neugebauer, Soifer, & Miley (1985), Iwasawa (1999), Vignad etal. (1999), and Sanders (2000). Although five of the X-ray detected submm sources are classified 
as AGNs, their submm emission is likely to be dominated by star formation. 


tor of the star-formation rate (e.g., Blain etal. 1999). However, 
five of the ten submm sources investigated here clearly host ob¬ 
scured AGNs. While the AGNs obviously dominate the X-ray 
emission in these sources, this does not necessarily imply that 
the submm emission is also powered by the AGNs. Unfortu¬ 
nately, the broad-band submm emission itself cannot provide an 
unambiguous diagnostic of the nature of the dominant physical 
mechanism since it is likely to be reprocessed thermal emission 
produced via the absorption of the primary emission by dust 
grains. 12 

CO observations can provide some insight into the origin of 
the submm emission by measuring the molecular gas mass and 
hence constrain the expected dust mass and emission from star 
formation (e.g., Young & Scoville 1991). For instance, CO 
observations have been used to show that star-formation activ¬ 
ity may account for « 50% of the submm emission from the 
z = 2.80 BALQSO SMM J02399-0136 (e.g., Frayer et al. 1998) 
and perhaps all of the submm emission from the z = 2.56 galaxy 
SMM J14011+0252 (e.g., Frayer etal. 1999; Ivison etal. 2001). 
However, in the absence of CO observations we must adopt an 


indirect strategy in estimating the origin of the submm emission 
from our submm sources. 

4.1. Constraints from the Submm-to-X-ray Spectral Slopes 

The relative importance of AGN and star-formation activ¬ 
ity in our submm sources can be estimated by comparing their 
submm-to-X-ray spectral slopes (asx) to those of well-studied 
nearby galaxies. We calculate asx as 

asx = - [log(#^) - 0.18] x 0.163, (3) 

7850^111 

where / 2 k e v is the observed flux density at 2 keV 
(keV cm' 2 s 1 keV 1 ), calculated from the full-band flux us¬ 
ing the estimated photon index, and /sso^m is the observed flux 
density at 850 pm (mJy). Larger values of asx indicate stronger 
submm emission relative to the X-ray emission. 

Following Fabian et al. (2000) we have compared our sources 
to the QSO 3C273, the Compton-thick AGN NGC 6240, and 
the archetypal ultra-luminous starburst galaxy Arp 220. While 
3C273 is clearly a powerful unobscured QSO, NGC 6240 con- 


12 Further insight into the origin of the submm emission may come from determining its spatial distribution (i.e., AGN emission is typically confined to the central 
galactic regions, while star formation can be more widespread). Such observations should be possible with the advent of the Atacama Large Millimeter Array (ALMA; 
see http://www.alma.nrao.edu/) in Chile and the Sub-Millimeter Array (SMA; see http://sma-www.harvard.edu/) in Flawaii. 
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tains components of obscured AGN and star-formation activ¬ 
ity. By contrast, the emission from Arp 220 appears to be en¬ 
tirely dominated by star formation. If an AGN is present in 
Arp 220, then to remain undetected at X-ray energies it must 
be either of low luminosity, or have extreme absorption fi.e., 
Ah ~ 10 25 —10 26 cm~ 2 ) and a small component of scattered 
emission (e.g., Iwasawa et al. 2001,2002; Clements et al. 2002). 
Since NGC 6240 has a considerably flatter X-ray spectral slope 
than 3C273 and Arp 220, the analysis of the submm-to-X-ray 
spectral slopes requires knowledge of the X-ray photon index 
for each source. For instance, at a given submm luminosity, the 
2 keV flux density of a T = 0.0 source will be m 5 times lower 
than that of a T = 2.0 source. Therefore, when interpreting the 
data for an unobscured QSO, the comparison should be made to 
3C273, while when interpreting the data for an obscured AGN, 
the comparison should be made to NGC 6240. 

The submm-to-X-ray spectral slope versus redshift of all the 
sources is shown in Figure 7. The AGN-classified sources have 
larger asx values than NGC 6240, suggesting that their AGNs 
are comparitively weaker at X-ray energies (with respect to 
their submm emission). Under the assumption that the obscura¬ 
tion in these sources is Compton thin (as expected for the major¬ 
ity of these sources; see §3.3), they are « 15-150 times less lu¬ 
minous at X-ray energies (for a given submm luminosity) than 
NGC 6240 (assuming a column density of Ah = 5 x 10 23 cm 2 ; 
i.e., the top dotted line in Figure 7). Taking all of the AGN- 
classified sources together, they are on average ss 70 times 
less luminous at X-ray energies (for a given submm luminos¬ 
ity) than NGC 6240. As can be seen from Figure 7, these re¬ 
sults are relatively insensitive to redshift so long as the AGN- 
classified sources lie at z > 1 • The only source where we have 
evidence that this may not be the case is the possible Compton- 
thick AGN 123622.6+621629 (see §3.3.1 and Figure 6). Under 
the assumption that the obscuration to this source is the same 
as that found for NGC 6240 (i.e., Compton-thick absorption 
with a 5% scattering fraction; Vignati etal. 1999), at z ~ 0.46 
it would be « 80 times less luminous at X-ray energies (for a 
given submm luminosity) than NGC 6240 (see Figure 7). 

We suggested in §3 that the X-ray emission from the submm 
sources without good X-ray constraints may be produced by 
star-formation activity in some cases. This hypothesis is fur¬ 
ther supported for the two sources with millimetric redshifts of 
Z <4 since their submm-to-X-ray spectral slopes are compara¬ 
ble to that expected from Arp 220 at z > 1 (see Figure 7). How¬ 
ever, we cannot distinguish between star-formation and AGN 
scenarios for the three z> 4 sources. 

4.2. Constraints from NGC 6240 

Although the individual constraints on the origin of the 
submm emission from the AGN-classified sources are weak, 
they all clearly have larger submm-to-X-ray spectral slopes 
than NGC 6240. Even though NGC 6240 is a nearby source 
(« 115 Mpc in our chosen cosmology), the origin of its submm 
emission is poorly known. For instance, the application of ra¬ 
diative transfer models shows that the far-IR-to-submm emis¬ 
sion can be almost completely explained by starburst activity 
(Rowan-Robinson & Crawford 1989), and the infrared (imag¬ 
ing and spectroscopy) characteristics of NGC 6240 suggest that 
star formation dominates the bolometric luminosity (Genzel 
etal. 1998). Furthermore, its nuclear optical spectrum shows 
the characteristics of a LINER (e.g., Veilleux et al. 1995). Con¬ 
versely, consideration of the energetics of the AGN derived 


from X-ray observations suggests that it may contribute « 50% 
of the far-IR-to-submm luminosity (e.g., Iwasawa & Comastri 
1998; Vignati et al. 1999; Lira et al. 2002). 

Assuming that the AGN completely dominates the submm 
emission in NGC 6240 (an extreme hypothesis) and the X-ray 
emission from the AGN scales linearly with submm luminos¬ 
ity, the contribution to the submm emission from AGN activ¬ 
ity in the AGN-classified sources would be 0.6-6.6% (as deter¬ 
mined from the « 15-150 times difference in X-ray luminosity 
for a given submm luminosity; see §4.1). Following the same 
assumptions, the average contribution to the submm emission 
from AGN activity would be « 1.4%. Hence, although there 
is considerable uncertainty in the AGN contribution at submm 
wavelengths, these analyses strongly suggest that the submm 
emission is dominated by star-formation activity. This conclu¬ 
sion is likely to be valid even if the z> 4 sources are also found 
to contain AGNs (see Figure 7). 

5. DISCUSSION 

It has been known for some time that AGNs comprise a frac¬ 
tion of the submm source population. Indeed, the first iden¬ 
tified source in a SCUBA blank-field survey was found to be 
an AGN (SMM J02399-0136; Ivison etal. 1998), and at least 
three of the 15 submm-detected sources in the SCUBA-lens 
survey are AGNs (Smail et al. 2002). In all these cases the ev¬ 
idence for AGN activity was found via optical spectroscopy 
of optically bright sources. The advantage of searching for 
AGNs with X-ray observations is that the limitations of optical 
spectroscopy are removed, allowing optically faint AGNs and 
AGNs without strong emission lines to be identified (e.g., Vi¬ 
gnati etal. 1999; Mushotzky etal. 2000; Alexander et al. 2001; 
Hornschemeier et al. 2001). However, obviously the lack of op¬ 
tical spectroscopy for the majority of our sources leads to con¬ 
siderable uncertainties in their source redshifts. In this final 
section, we discuss the dependence of our results on the as¬ 
sumed source redshifts, the fraction and properties of AGNs in 
the submm source population, and the AGN contribution to the 
submm emission in SCUBA sources. 

5.1. Dependence of the Results on Source Redshifts 

The redshifts for the majority of the submm sources are de¬ 
rived from the radio-to-submm spectral slope and are uncertain. 
As mentioned in §4.1 these uncertainties should not change 
our main result (i.e., that the average AGN contribution to the 
submm emission is of the order of « 1.4%) so long as the 
sources lie at z > 1. However, we note here that since the red¬ 
shift tracks for the template galaxies rise to lower values of asx 
at z < 1 (see Figure 7), the case for star formation dominating 
the submm emission is even stronger for sources at lower red¬ 
shifts (as may be the case for 123622.6+621629; see §3.3.1). 
Finally, we also note that although the implied column densi¬ 
ties and X-ray luminosities of the AGN-classified sources will 
be lower if they lie at z < 1, since T < 1.0 power-law emission 
is only likely to be produced by obscured AGN activity, their 
classification as AGNs is insensitive to redshift. 

5.2. AGNs in the Bright Submm Source Population 

In this study we have shown that five of the seven X-ray de¬ 
tected submm sources appear to be obscured AGNs, placing 
a lower limit on the source density of obscured AGNs with 
> 5 mJy of 260+}™ deg 2 . Under the same assump¬ 
tions as those used in determining the fraction of the bright 
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submm source population with X-ray counterparts (see §2.4.2), 
this suggests that £ 26% of the bright submm source popu¬ 
lation host AGNs (in reasonable agreement with model predic¬ 
tions; e.g., Almaini, Lawrence, & Boyle 1999; Fabian & Iwa- 
sawa 1999). 

Ivison et al. (2002) found a similar AGN fraction to that 
derived here using the moderately deep XMM-Newton ob¬ 
servations of the Lockman Hole East region (e.g., Hasinger 
etal. 2001). However, their submm-detected AGN source 
density is w 2 times smaller (i.e., 120+jo deg~ 2 versus 
260+[™ deg~ 2 ). 13 The difference in the submm-detected AGN 
source densities is mostly due to the higher sensitivity of our 
X-ray observations. Hence the similarity between the AGN 
fraction found by Ivison et al. (2002) and our lower limit may 
indicate that our estimated AGN fraction is conservative. 

The unabsorbed rest-frame full-band luminosities of the 
AGN-classified sources are Lx ~ 10 43 —10 44 erg s~ 1 . This 
range of X-ray luminosity is typical of luminous Seyfert galax¬ 
ies. Therefore, although the majority of the submm sources 
are likely to be extremely luminous infrared galaxies (e.g., Ivi¬ 
son et al. 1998; Smail et al. 2002), their AGNs are comparitively 
weak. For example, typical QSOs of similar infrared luminos¬ 
ity (e.g., Isaak et al. 2002; Page etal. 2001) are « 10-100 times 
more luminous at X-ray energies (e.g.. Page etal.2001; Vig- 
nali et al. 2001,2003). Under the assumption that bright submm 
sources are massive proto galaxies (e.g., Blain etal. 1999; Lilly 
etal. 1999), the luminosities of these AGNs are consistent with 
those expected for growing supermassive black holes (e.g., 
Archibald et al. 2002). 

5.3. AGN Contribution to the Submm Emission 

The current constraints suggest that the total contribution to 
the submm emission from AGN activity is unlikely to be more 
than « 1.4%; however, as noted in §2.2, there are likely to be 
further bright submm sources in our region that have not yet 
been detected by SCUBA. In addition, we have not considered 
the possible contribution from optically bright submm-detected 
QSOs. Although optically bright submm-detected QSOs may 
comprise only « 5-10% of the bright submm source population 
(see footnote 1), their submm emission could have a large AGN 
component (perhaps of the order of « 50%; Frayer etal. 1998; 
Bautz etal.2000). Hence, although optically bright submm- 
detected QSOs may not be numerically significant, they may 
provide the bulk of the AGN emission at submm wavelengths. 

6. CONCLUSIONS 

We have used the 2 Ms Chandra exposure of the CDF-N 
to constrain the X-ray properties of bright SCUBA sources 
(/ 850 /jm > 5 mJy; S/N> 4). In this study we have focused 
on the X-ray spectral properties of the X-ray detected submm 
sources to determine whether the AGN-classified sources are 
Compton thick or Compton thin. We have used these results 
to constrain the contribution to the submm emission from AGN 
activity. Our main results are the following: 

• Seven of the ten bright submm sources are detected 
with X-ray emission. The corresponding source den¬ 
sity of bright submm sources with X-ray counterparts 
( 36 O+J 3 Q deg -2 ) suggests that ^ 36% of bright submm 


sources have X-ray counterparts at this X-ray depth; 
however we note that this fraction is somewhat uncer¬ 
tain since this region does not have complete SCUBA 
coverage down to = 5 mJy. See §2. 

• Five of the X-ray detected submm sources have flat X- 
ray spectral slopes (r < 1.0) and luminous X-ray emis¬ 
sion (L x « 10 43 - 10 44 erg s 1 ), suggesting obscured 
AGN activity with Ah A, 10 23 cm 4 . One source is 
possibly a Compton-thick AGN since it has an ex¬ 
tremely flat X-ray spectral slope (r w -0.7) and shows 
possible evidence (« 2.3 o) for a redshifted 6.4 keV 
iron Ka emission line; however, the inferred redshift 
(z = 0.46+QQ2) is in considerable disagreement with the 
millimetric redshift (z, = 2.4+Qg). See §3. 

• A comparison of the five AGN-classified sources to the 
well-studied heavily obscured AGN NGC 6240 sug¬ 
gests that the AGNs in these sources contribute, on av¬ 
erage, a negligible fraction (i.e., « 1.4%) of the submm 
emission. Hence, the submm emission from these 
sources appears to be dominated by star formation. This 
result is relatively insensitive to redshift. See §4. 

• The X-ray constraints are weak for the other five submm 
sources. We find that the X-ray properties for the two 
sources with millimetric redshifts of z < 4 are consis¬ 
tent with those expected from luminous star-formation 
activity (Lx « 10 42 -10 43 erg s _1 ). The submm-to-X-ray 
spectral slopes of these sources are also similar to those 
expected from Arp 220 (the archetypal ultra-luminous 
dusty starburst galaxy) at z > 1 ; however, we cannot rule 
out the possibility that low-luminosity AGNs are present 
in these sources. The three z> 4 sources could be either 
AGNs or starbust galaxies. See §3 and §4. 

• The fraction of the bright submm source population with 
AGN activity (i.e., £ 26%) is in reasonable agreement 
with model predictions. The comparitively low X-ray 
luminosities of the AGN-classified sources are more 
consistent with Seyfert galaxies than QSOs. We suggest 
that optically bright submm-detected QSOs (i.e., lumi¬ 
nous type 1 AGNs), although possibly not numerically 
significant, may provide the bulk of the AGN emission 
at submm wavelengths. See §5. 

These results show that the low X-ray detection rate of bright 
submm sources by moderately deep X-ray surveys is proba¬ 
bly due to the relatively low luminosities of the AGNs in these 
sources rather than Compton-thick absorption. The current con¬ 
straints therefore suggest that the total contribution from AGNs 
to the submm emission in the bright submm source population 
is negligible. Hence, the submm emission of bright submm 
sources appears to be predominantly powered by star-formation 
activity, and therefore can be used to determine star-formation 
rates. 


13 Bautz etal. (2000) also found a similar AGN fraction from combining their two Chandra detected submm sources with the Chandra undetected submm sources in 
the HDF-N (e.g., Homschemeier et al. 2000), although their AGN fraction may be biased by the inclusion of a comparitively rare submm-detected QSO (the BALQSO 
SMM J02399-0136; Vemet & Cimatti 2001). 
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Table 1 

Basic properties of the bright submm sources 


Coordinates 

Radio 

X-ray 



Radio 




Submm 

Derived properties 
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^2000 3 

detected? 
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a Source coordinates. Taken from, in order of preference, the radio position (if radio detected), the X-ray position (if X-ray detected), or the submm source position. 

b Radio flux density at 1.4 GHz in units of ^uJy and radio spectral index (a). The radio flux density upper limits correspond to the 5 a survey limit given in Richards 
(2000). The radio spectral index is defined from the 8.5 GHz and 1.4 GHz flux densities as F v ©< v -a . Taken from Richards (2000) and Dunlop etal. (2002). 

c Source classification based on radio properties for sources detected in Richards (1999). “S” indicates star-forming galaxy, and “U” indicates that good constraints could 
not be placed. 

d /-band Vega-based magnitude and I — K color. The calculated I — K color is determined from I — (HK' — 0.3) following Barger et al. (1999b) for all sources except 
123652.0+621225. The /-band magnitude and I — K color for 123652.0+621225 are taken from Dunlop et al. (2002). 

e Submm flux density at 850 /um in units of mJy. 

Reference for submm data (HU=Hughes et al. 1998; B00=Barger et al. 2000; B02=Barger et al. 2002; BY=Borys et al. 2002). 

g Source redshift, either spectroscopic (Cohen et al. 2000; Hornschemeier et al. 2001), based on the source multi-wavelength SED (Dunlop et al. 2002), or millimetric 
(Carilli & Yun 2000); only 123629.1+621045 has a spectroscopic redshift (z = 1.013) and only 123652.0+621225 has an SED-based redshift ( z = 4. 11q. g). The uncertainties 
in the SED and millimetric redshifts correspond to the la level. We note that the source redshift of 123622.6+621629 may be z ~ 0.46, based on X-ray spectral analysis 
(see §3.3.1). 

h Rest-frame 1.4 GHz radio luminosity density in units of 10 24 W Hz -1 determined following equation 2. The uncertainties on the calculated luminosity densities are 
determined using the upper and lower redshift bounds. We are unable to place luminosity constraints for the three radio-undetected sources as they all have redshift lower 
limits. 

'This source may have been detected at low-significance in the 1 Ms CDF-N study of VROs (Alexander et al. 2002b); see §2.4. 
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Table 2 

X-RAY PROPERTIES OF THE BRIGHT SUBMM SOURCES 


Coordinates 


X-ray counts 


Band 

Effective 


X-ray flux 


Derived properties 

X-ray 

0^2000 3 

§ 2000 a 

FB b 

SB b 


Ratio c 

pi 

FB C 

SB C 

HB e 

<*sx f 

L x t 

ID" 

12 36 16.15 

+62 15 13.7 

i3o.+;« 

68 .8i; 0 „ 2 

57.4i;°3 8 

0.851°;?° 

0 Q+ 0-2 
uy - 0.2 

1.01 

0.16 

0.80 

1 to +002 

8.3+“/ 

AGN 

12 36 18.33 

+62 15 50.5* 

29.4 

<17.5 

<17.2 

- 

1.4 

0.17 

<0.04 

<0.19 

1 41 + 0.03 
1 '^ 1 -0.03 

0.41“ 

U 

12 36 20.3 

+62 17 01 

<15.5 

<10.4 

<14.1 



<0.10 

<0.03 

<0.19 

>1.49 


U 

12 36 21.1 

+62 12 50 

<11.1 

<5.3 

<10.3 



<0.07 

<0.01 

<0.14 

>1.50 

- 

u 

12 36 22.65 

+62 16 29.7 

63.41“;? 

<13.2 

56.9l“ 4 7 

>4.39 

<-0.6 

1.03 

<0.03 

1.15 

1 4 _ | +0.03 
i- ^-0.03 

13.2lg 3 

AGN 

12 36 29.13 

+62 10 45.8 

198.4TJH 

66.31“,° 

134.71&J 

2 05 +0 - 38 

0.33 

p 

1 "+ 

p p 

2.32 

0.16 

2.26 

1 28 +003 
0.03 

3.1 

AGN 

12 36 46.05 

+62 14 48.7 

13.8+JS 

12 .6l|| 

<7.5 

<0.60 

1.4 

0.08 

0.03 

<0.09 

1 4Q+0-06 

0.41“ 

U 

12 36 52.06 

+62 12 25.7 

<9.8 

<10.8 

<4.0 



<0.06 

<0.03 

<0.05 

>1.48 

<1.1 

u 

12 37 07.21 

+62 14 08.1 

«4 4+H - 4 
°^- 10.2 

32.6ip 

52.31” 

1 62 +0 ’ 48 
I-OZ—0 39 

0 T +0 - 3 
u>: 5 - 0.2 

0.98 

0.09 

0.92 

1 ->4+0.03 
1 " >t -0.03 

1 c q+24.5 

AGN 

12 37 12.07 

+62 12 11.3 

30.71” 

10.512? 

20.91” 

2 00+ 145 
Z.UU—0 88 

0 l +0-5 

ul -0.5 

0.39 

0.03 

0.38 

1 A A +0-02 

1 ‘^- 0.02 

>11.3 

AGN 


a Source coordinates. Taken from Table 1. 

b Source counts and errors or upper limit. Determined following the procedure given in Brandt et al. (2001) for sources detected with WAVDETECT false-positive 
probability threshold of 10 -7 and from WAVDETECT for sources detected with WAVDETECT false-positive probability threshold of 10 -5 . “FB” indicates full band, “SB” 
indicates soft band, and “HB” indicates hard band. Upper limits are calculated following §3.2.1 of Brandt et al. (2001). 

c Ratio of the count rates in the 2.0-8.0 keV and 0.5-2.0 keV bands. The errors were calculated following the “numerical method” described in §1.7.3 of Lyons (1991). 

d Effective photon index for the 0.5-8.0 keV band, calculated from the band ratio. The photon indices for the X-ray sources detected with a threshold of 10 -5 and those 
with a low number of counts have been set to T = 1.4, a value representative of the X-ray background spectral slope; see Brandt et al. (2001). The photon index is related 
to the energy index by a = T-l where F v <x v -a . 

e Fluxes are in units of 10 -15 erg cm -2 s -1 . These fluxes have been calculated following the method described in Brandt et al. (2001). They have not been corrected for 
Galactic absorption. “FB” indicates full band, “SB” indicates soft band, and “HB” indicates hard band. 

f Submm-to-X-ray spectral index (CCsx) calculated following equation 3. 

g Rest-frame full-band luminosity in units of 10 43 erg s -1 calculated following equation 1. The uncertainties in the calculated luminosities are determined using the upper 
and lower redshift bounds. For the AGNs the unabsorbed luminosity is given (corrected assuming the absorption is Compton thin and the underlying unabsorbed emission 
is a r = 2.0 power law), while for the other sources T = 2.0 is assumed. We are unable to place luminosity constraints for the two X-ray-undetected sources with redshift 
lower limits. 

h Source classification based on X-ray properties; see §3.2. “AGN” indicates an AGN and “U” indicates that good X-ray spectral constraints could not be placed. 

‘X-ray counterpart detected with WAVDETECT with false-positive probability threshold of 10 -5 . 
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Table 3 

X-RAY SPECTRAL FITS FOR THE AGN-CLASSIFIED SOURCES 


Object name 

PL 

PL+ABS 


r a 

r b 

Nft (cm 2 ) 

123616.1+621513 

1 o +a3 

lu -0.3 

o 9+U.S 
^ -0.7 

0-8^4 

123622.6+621629 

1 

O 

Li 

i + 
o c 

4^ 4^ 

0 1+ 10 
u - 1 -0.8 


123629.1+621045 

0 3+ 0 - 2 
U - J -0.2 

i o+0.6 

1 * —0.5 

i o +0 - 6 

i ' U -0.5 

123707.2+621408 

0 5+ 0 - 3 
u - J -0.3 

i-6±8:® 

1 l+°- 9 

123712.0+621211 

o o+ 0 - 5 
U ' U -0.6 

0.5±?;1 

<0.8 


NOTE. — The uncertainties refer to the 90% confidence level (for one interesting parameter). 
a Best-fit photon index for Galactic absorption model. 
b Best-fit photon index for absorbed power-law model. 

c Best-fit column density in units of 10 22 cm -2 at z — 0; the column density at the source redshift (A/h, z ) is related to the best-fit column density at z — 0 by 

Wh,z«(1+z) 2 ' 6 /Vh. 



